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We show that correlations inconsistent with any locally causal description can be a generic feature 
of measurements on entangled quantum states. Specifically, spatially-separated parties who perform 
local measurements on a maximally-entangled state using randomly chosen measurement bases can, 
with significant probability, generate nonclassical correlations that violate a Bell inequality. For n 
parties using a Greenberger-Horne-Zeilinger state, this probability of violation rapidly tends to unity 
as the number of parties increases. We also show that, even with both a randomly chosen two-qubit 
pure state and randomly chosen measurement bases, a violation can be found about 10% of the 
time. Amongst other applications, our work provides a feasible alternative for the demonstration of 
Bell inequality violation without a shared reference frame. 
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One of the most remarkable features of quantum me- 
chanics is that distant measurements can exhibit corre- 
lations that are inconsistent with any locally causal de- 
scription (LCD) These quantum correlations are sig- 
natures of entanglement ||, and serve as a resource || 
for a range of classically-impossible information process- 
ing tasks such as quantum key distribution Q, tclcpor- 
tation JBj and reduced communication complexity j6| . 

Not all quantum states are useful for quantum in- 
formation processing tasks. Even for entangled states, 
nonclassical correlations are not an immediate conse- 
quence of the entanglement present in such systems (see 
Rcf. JtJ and references therein). They also depend cru- 
cially on the choice of measurements to which they are 
subjected, and specifically how these measurements are 
correlated between parties. For example, in the stan- 
dard scenarios to violate a Bcll-Clauser-Hornc-Shimony- 
Holt (Bcll-CHSH) inequality |) or more generally an 
n-party inequality using a Greenberger-Horne-Zeilinger 
(GHZ) state, a maximum violation is achieved by all par- 
ties using specific measurements in the common x — y 
plane ||. In an experiment, this requires some care [ [To| ; 
for two parties using polarization-encoded single-photon 
pairs transmitted through long bircfringcnt optical fibers, 
the random rotation of the polarization of each photon is 
first compensated through some method of alignment in 
order to set this frame jnj . With larger numbers of par- 
ties, greater violations can be achieved but more parties 
require an increasing complexity of alignments. 

Several approaches to circumvent this problem have 
been studied. Encoded entangled states that are in- 
variant with respect to some collective unitary opera- 
tions |l^, [l3| can be used, but generally require much 
more complicated state preparation as well as joint mea- 
surements on multiple spins. Alternatively, an alignment 
of frames can be performed through the coherent ex- 
change of quantum systems O, |l4| or by supplement- 



ing each system with a small quantum reference frame of 
bounded size JD|. Such solutions are resource- intensive 
given that this alignment consumes many of these quan- 
tum resources that could otherwise be used to generate 
useful correlations. 

Here, we prove that such resource- intensive approaches 
are unnecessary, and that nonclassical correlations oc- 
cur with high probability even without any alignment of 
local measurements. Specifically, we investigate the de- 
tection of nonlocal correlations wherein each party uses 
randomly chosen measurement bases. Ruling out LCD 
in this scenario cannot be achieved dctcrministically for 
a general entangled quantum state [[l3|. However, we 
demonstrate that the probability of observing correla- 
tions inconsistent with LCD using randomly chosen mea- 
surement bases can be remarkably high. Although one 
might naively expect that increasing the number of par- 
ties participating in the test (and the corresponding com- 
plexity of alignments) would only complicate matters, we 
show that the probability of demonstrating a violation 
rapidly approaches unity as the number of parties in- 
creases. Finally, we demonstrate the possibility of two 
parties violating a Bell inequality using both randomly 
chosen measurements and a randomly chosen pure state. 

We consider a typical scenario of an n-party Bell- 
experiment using spin-1/2 particles, wherein a verifier 
prepares N 3> 1 copies of the n-partite GHZ state 

|GHZ) n = ^(|0) 1 |0) 2 ...|0) n + |l) 1 |l) 2 ...|l) n ) ! (1) 

and distributes them to n parties. Each particle is sub- 
jected to a local measurements with two possible out- 
comes, ±1. The measurement for each particle is cho- 
sen randomly, with equal probability, from a set of two 
measurement bases, defined for the kth party by a pair 
of spatial directions fil^' (0, </>), with Sk = 1,2 labeling 
the bases. For spins measured via a Stern-Gerlach de- 
vice, these directions correspond to the orientation of the 
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magnetic field in the measurement apparatus. However, 
an equivalent picture in terms of spatial directions can 
be used for any quantum systems described by a two- 
dimensional Hilbert space (eg., the polarization encoding 
of single photons) via the Bloch sphere. After the com- 
pletion of all measurements, the n parties return to the 
verifier a list containing only the label Sk of which mea- 
surement they performed for each particle and the value 
±1 of the outcome they obtained for that measurement. 

A relevant Bell inequality for this scenario is the n- 
partitc Mermin-Ardchali-Bclinskii-Klyshko (MABK) in- 
equality [pjf- 
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Here, the n-partite correlation function E(si, . . . , s n ) is 
defined to be the expectation value of the product of 
measurement outcomes (±1) given that the fcth party 
measured their spin in the basis When n = 2, the 
MABK inequality is equivalent to the familiar Bcll-CHSH 
inequality 5 CH sH = \E{1, l)+E(\, 2)+E(2, 1)-E(2, 2)| < 

In a standard Bell experiment, the maximal violations 
for a given Bell inequality are achieved by choosing the 
local orientations and O^' to be perpendicular M, 
i.e., ■ Q 2 = 0, an< i by fixing the relative orientations 
of these measurements with the other experimentalists 
(and the verifier who prepared the systems). The latter 
necessitates the use of a shared reference frame; in what 
follows, we will sec what can be achieved without a shared 
reference frame. 

We will first consider the completely unconstrained 
case, wherein each party k chooses both members of their 
pair of directions for sj. = 1,2 independently and 
uniformly from the set of all possible directions. We re- 
fer to this case as random isotropic measurements (RIM) . 
In such a scenario, we determine the probability ^Liabk^ 
that the measurement statistics received by the verifier 
will give rise to a violation of the MABK inequality of 
Eq- (PI), given by 
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where dili^ / 4ir is the Haar measure associated with mea- 
surement direction fii* 1 and /" ABK ({^}£ = i) is a func- 
tion that returns 1 if orientations {Oi^}^ =1 give rise to 
a correlation (i.e., measurement statistics), that violates 



the MABK inequality and otherwise. These integrals 
simply pick out those measurement directions that will 
give a violation, normalizing them against the total vol- 
ume of possible measurement directions. 

For the detection of nonclassical correlation, we can 
allow the verifier to make use of the freedom in rela- 
beling all measurement settings and/or outcomes and/or 
parties jlT], [H]]. To appreciate the importance of this 
fact, consider as an example the following choice of 
local measurement for the n , = 2 case: {cr Zl cr x } for 
Alice and (a x + a z ) , -±= {a x - a z )} for Bob. Sup- 
pose that Alice has labeled her measurements such that 
a z <-> si=l and a x <-> si=2, while Bob has labeled 
his measurements such that =j= (o~ x + cr z ) <-> S2=l and 

^ (o~ x — o~ z ) f-> S2=2, then the verifier will find that as 
it is, <S C hsh vanishes, and hence does not lead to a Bcll- 
CHSH inequality violation. However, if the verifier is 
also given the freedom to relabel the measurement set- 
tings, they could use, instead, the following labeling for 
Alice's measurements: a x <H> Si=l and a z O si=2. As 
opposed to the original labeling, this new labeling un- 
veils measurement statistics that violate the Bcll-CHSH 
inequality maximally. We say that two Bell inequalities 
are equivalent if they can be obtained from each other 
by relabeling the measurement settings and/or outcomes 
and/or parties |l7], In calculating the probability of 
observing a violation, we clearly want to include all such 
equivalent inequalities. 

We numerically calculate Pmabk ^ ; using the relabeling 
strategy to test a given set of experimental data against 
all 2 n+1 equivalent MABK inequalities @; see Fig. |. 
The special case of n = 2 can be solved analytically (see 
supplementary materials attached below for details) and 
is found to be 2(ir — 3) ss 28.3% (see also Rcf. for con- 
trasting results). For n > 2, the probability is less than 
this amount but increases monotonically for n > 3 up to 
the limit of our analysis, n = 15. Note, however, it ap- 
pears that this probability is asymptotically approaching 
a value that is less than unity. 

The case of n = 2 appears somewhat anomalous. 
The reason is that, for n = 2, the equivalence class of 
Bell-CHSH inequalities obtained by relabeling 6>mabk ( or 
Scusu) is both necessary and sufficient for determining if 
a given set of experimental statistics generated by per- 
forming two binary-outcome measurements per site ad- 
mits a LCD ^2j. In other words, in the space of mea- 
surement statistics generated by performing such mea- 
surements, all nonclassical correlations can be detected 
by this equivalence class of inequalities. 

For n > 2, there are (i) Bell inequalities with n-partitc 
correlation functions that are not of the MABK form 
(^) and (ii) Bell inequalities involving less-than-n-partite 
correlation functions. To determine if a given correla- 
tion for n > 2 is incompatible with any LCD, we need 
to test the measurement statistics against the complete 
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FIG. 1: (Color online) Probability of finding a nonclassical 
n-partite GHZ correlation. The two sets of data at the bot- 
tom and top are, respectively, for correlation generated from 
random isotropic measurements (RIM) and random orthogo- 
nal measurements (ROM). The markers are for the scenarios 
when the verifier has access to (1) the 2 n+1 MABK inequal- 
ities (A for RIM, x for ROM), and (2) the 2 2 " WWZB in- 
equalities (V for RIM, + for ROM). 
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74.6899% 
5000 


94.2380% 
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99.5926% 
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41.2982% 
5000 


96.2073% 
5000 


99.9757% 
500 


99.9999% 
25 


100.00% 
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TABLE I: Probability of finding a nonclassical correlation 
from the n-partite GHZ state for the scenario where each 
party is allowed to perform binary projective measurements 
in two randomly chosen measurement bases. The number of 
data points (Run) used to compute the probabilities for RIM 
(ROM) is included in the third (fifth) row of the table. 



set of Bell inequalities that are relevant to the particu- 
lar experiment. For a small number of parties, this can 
be achieved efficiently using linear programming (sec ]23[ ] 
and the supplementary materials attached below) or by 
checking the correlation directly against the complete set 
of relevant Bell inequalities. Our results arc summarized 
in Table § 

These results clearly indicate that the probability of 
a violation rapidly approaches unity with increasing n. 
The input to the linear program, however, scales expo- 
nentially with n (see supplementary materials attached 
below), thus making it intractable to determine the prob- 
ability of violation reliably for larger values of n. To ob- 
tain some insight into the behavior, we can restrict our 
attention to the extensive set of 2 2 n-partite correla- 
tion inequalities discovered independently by Werner & 
Wolf and Zukowski & Brukncr (l| (WWZB). While 
violation of these inequalities is still sufficient to rule out 
LCD, the converse is generally not true. Thus, using this 



set of inequalities yields a lower bound on the probability 
that a given correlation is nonclassical. 

To test if a given correlation satisfies all of the WWZB 
inequalities is equivalent to testing if it satisfies the fol- 
lowing nonlinear inequality p7|, |l8| 



E E 



= ±1 si, 



=ij=i 



-\e( S1) 



< 2 r 



(5) 



We have numerically computed the probability of finding 
a randomly generated correlation from RIM to violate 
inequality (Eh for n < 15. As can be seen in Fig. [I], 
this probability generally increases with n and is above 
50% for n = 15. However, it is inconclusive whether this 
probability asymptotes to unity or not for large n. 

Finally, we consider another mechanism by which par- 
ties who do not share a reference frame can increase the 
probability of observing a violation. So far we have con- 
sidered each of the two local measurement bases for each 
party to be chosen independently and isotropically. How- 
ever, because the optimal violations are typically found 
when the pair of local measurement directions are or- 
thogonal |§, 

we can consider a random selection of mea- 
surement bases under the constraint that each local pair 
are orthogonal. We have repeated the above calculations 
with this constraint (i.e., that • fi^ = 0); we refer 
to this distribution as random orthogonal measurements 
(ROM). The results reveal a significant increase in the 
probability of a violation, as shown in Fig. [l] and Table Q. 
For larger values of n, our results for the WWZB inequali- 
ties also strongly suggest that the probability asymptotes 
to unity for large n. 

Our results so far have all made use of a fixed 
maximally-entangled pure state (Q) (although clearly the 
local bases in which this state is defined are irrelevant 
for our result). We can also consider how the probability 
of observing a violation of LCD varies with the degree 
of entanglement of the distributed quantum state. For 
n = 2, a pure quantum state can always be written as 
|*) = cos#|0)i|0)2 + sill 6>| 1) 1 1 1) 2 for some local bases. 
Numerically, we can compute the probability of viola- 
tion as a function of entanglement (for details, see sup- 
plementary materials attached below). With this result, 
we can determine the probability of observing a viola- 
tion of LCD using both randomly chosen measurement 
bases and a randomly chosen state. We sample two-qubit 
pure states randomly and uniformly from C 2 <8>C 2 accord- 
ing to the SU(4) Haar measure. Given that our results 
are independent of local choice of bases, we can repro- 
duce this distribution by sampling Schmidt coefficients 
as for single-qubit mixed states from the uniform Bloch 
ball p3]. The relevant measure is the Hilbert- Schmidt 
measure on a single qubit, i.e., P HS (r) = 24r 2 , where 
r = 
Ref. 

chosen pure two-qubit state can then be determined by 



^|2£ i s the length of the Bloch vector (pp. 354, 
p4|). The probability of violation for a randomly 
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numerical integration to be, respectively, about 5.3% for 
RIM and 10.1% for ROM. 

Conclusion and other directions - We introduced the 
idea of performing a Bell experiment using randomly- 
chosen measurement bases, with the aim of demonstrat- 
ing nonclassical correlations in the absence of a shared 
reference frame. We applied this idea to the n-partitc 
GHZ state and showed that the probability of finding 
a nonclassical correlation rapidly tends to unity as the 
number of parties increases. Naively, one might have ex- 
pected that this chance would diminish with increasing 
n, due to the increasing chance of misalignment of the 
experimenters' measurement bases (or, alternatively, by 
considering the fragility of the GHZ state to dephasing). 
Our results clearly show otherwise. We have thus showed 
that without a shared reference frame, a Bell inequality 
violation can still be demonstrated reliably without re- 
sorting to complicated state preparation or consumption 
of expensive quantum resources. This, of course, sig- 
nificantly reduces the technical requirements for experi- 
mentally violating a Bell inequality, performing quantum 
key distribution based on such violations, or establishing 
large-scale quantum networks. 

Our work also represents the first systematic study of 
the set of correlation derivable from a quantum resource 
such as the n-partitc GHZ state. Specifically, our re- 
sults indicate that the correlations obtained by perform- 
ing projective measurements on this quantum resource 
are almost ubiquitously nonclassical (this can be rigor- 
ously quantified by considering a volume measure in the 
space of GHZ-attainable-corrclations induced by RIM 
and ROM). Given that measurements and the resulting 
correlations play an important role in many quantum 
information processing tasks, the tools and ideas that 
we have introduced here, suitably generalized, may shed 
some light on the origin of the power of other quantum re- 
sources such as the highly entangled cluster state. In par- 
ticular, they may add new insights to the computational 



power of correlations recently discussed in Ref. 25 



The results presented here motivate several additional 
research directions. First, one may attempt to determine 
a new intuition for our result of increasing probability of 
violation with the number of parties n, by relating this 
probability to the average magnitude of the MABK vi- 
olation, which also increases with n. Second, for many 
quantum information processing tasks, violation of the 
inequalities considered in the present work may not be 
sufficient, as other types of multipartite entanglement 
(perhaps different from that of the GHZ state) may be 
required that are not quantified by these inequalities. 
To this end, it will be interesting to investigate the effi- 
ciency of random measurements against inequalities that 
detect other, more general types of multipartite entan- 
glement [^6| . Finally, for real- world applications of these 
results, a critical issue is to determine how robust they 
are against decoherence. 
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SUPPLEMENTARY MATERIALS 
MABK-violation without classical post-processing 

In the main text, we have presented the probabil- 
ity of observing a violation by considering all equiv- 
alent Mermin-Ardehali-Belinskii-Klyshko (MABK) in- 
equalities. If the classical post-processing to optimize 
over the measurement labelings is not carried out (i.e., we 
calculate the probability to violate a particular inequality 
rather than an equivalent set), then wg find Pmabk 

; de- 
creases steadily with n; see Fig. |2[ However, in this case, 
among those choices of bases that give rise to an MABK- 
violation, we find that their average violation increases 
with n when normalized by the classical threshold 2™~2 
but decreases when normalized by the maximal possible 
quantum violation 2t _1 [|[(| [l7| . See also Ref. p7| in 
this regard. 
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FIG. 2: (Color online) Probability of finding a randomly 
generation correlation from |GHZ) n to violate a single MABK 
inequality. 



Analytic derivation of probability of violation for 

n = 2 



For the specific case of n = 2, the probability of violat- 
ing the (2-party) MABK inequality, i.e., the Bell-CHSH 
inequality § § 

Schsh = \E{1, 1) + E(l, 2) + E{2, 1) - E(2, 2)| < 2, (6) 
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can be computed analytically. To this end, it is expedient 
to consider the two-qubit singlet state 

l*-> = ^(|0>i|l)2-|l)i|0) a ), (7) 

instead of the 2-party GHZ state: 



|GHZ> 2 = ^(|0>i|0) 2 + |l) 1 |l> 2 ). 



(8) 



Clearly, these states are local unitarily equivalent, and 
hence their probability of violating the Bell-CHSH in- 
equality via random isotropic measurements (ROM) is 
identical. However, \^~) is also rotationally invariant 
(i.e., U^)U\^f~) = \^~) for arbitrary qubit unitary trans- 
formation U), which allows us to simplify our analysis 
considerably. 

Let us begin by writing the measurement directions 
for Alice and Bob as a si and b s . 2 . From some simple 
calculation, one can show that the correlation function, 
i.e., the average of the product of measurement outcomes 
for Alice measuring a Sl ■ a and Bob measuring b Ss ■ 3 is 



(9) 



E [a Sl ,b S2 = -a Sl ■ b S2 , 



where a = (a x , <jy, <j z ) is the vector of Pauli matrices. 
Using this in Eq. (g), we are thus interested to know how 
often the inequality 



tii ■ b\ + tii ■ b 2 + 3 2 ■ b\ — ti 2 ■ b 2 > 2. (10) 



holds when the unit vectors ai, ti 2 , &i and 62 are cho- 
sen independently, randomly and isotropically. We can 
simplify this by defining the two orthogonal vectors, 
fj_ = bi — 62 and r = bi + b 2 , in terms of which in- 
equality (not) holds if and only if 



\tii ■ r + ti 2 ■ r±\ > 2. 



(11) 



Note that as well as being orthogonal, the vectors r± and 
r have lengths satisfying, 



\r\ = ^2(1 + 61 • b 2 ) = ^2(1 + as), 
|r±| = V // 2(l-&i-fo 2 ) = v / 2(l-x), (12) 

where x — bi ■ b 2 . 

By rewriting r in terms of its unit vector f and mag- 
nitude \f\, i.e., r = f\r\ = fy2(l + x), and similarly for 
rj_, inequality ([ll]) becomes 

\tii ■ f\f\ + x + ti2 ■ rj_\/l — x\ > V2. (13) 

Now the rotational invariance of dot products ensures 
that in determining the fraction of measurement direc- 
tions that violate inequality ([13]), the actual direction of 
f and fj_ is irrelevant. With some thought, one can see 



that to determine this fraction, and hence the probabil- 
ity of violating the Bell-CHSH inequality (0) via random 
isotropic measurements, it suffices to sample x and the 
dot products a\ • r, a% • fj_ uniformly from the interval 
[—1, 1]. Writing j3 = ai ■ f and a = ti 2 ■ f±, we there- 
fore wish to find the fraction of values of a, (3, x € [—1, 1] 
satisfying, 



laVl -x + /Vl + x\ > V2. 



(14) 



This problem has a useful geometrical interpretation. 
The set of all points a,(3,x e [—1,1] define a cube, <£, 
which encloses a volume 2 3 and has faces located at 
a, f3, x = ±1. To find the fraction of Alice and Bob's mea- 
surement directions that would violate the Bell-CHSH 
inequality we need to find what fraction of the volume 
of €. contains points (a,/3,x) satisfying the constraint of 
Eq. ©. 

By separately considering the cases where a\/l — x + 
/3\/I + x > and ay/1 — x + f3\/l + x < 0, the regions of 
<£ containing points satisfying Eq. (|l4|) are given by, 



and, 



a > 



a < 



V2-/VT 



V2 + I3VT 

y/T^X~ 



(15) 



(16) 



If we take a cross-section through £ corresponding to 
a particular fixed value of x, then these two expressions 
define triangular regions at the corners of the cube, as 
is shown in Fig. ||. From some straightforward compu- 
tation, it can be shown that the total area of the two 
triangles for a given x is 



vT 



,•2 



(17) 



Integrating this expression from x = — 1 to 1, and divid- 
ing the result by 2 3 (the total volume of £), we therefore 
find that the fraction of Alice and Bob's measurement 
directions that will violate the Bell-CHSH inequality is: 



7T-3 



7.0796%. 



(18) 



Notice that the above expression is the probability of 
violating inequality (||) without taking into account of 
any possible relabeling of measurement settings and/or 
outcomes. To take this into account, note that the four 
equivalent Bell-CHSH inequalities are as follows: 

\E(1, 1) + E(l, 2) - E(2, 1) - E(2, 2)| < 2, (19) 

\E(1, 1) + E(l, 2) - E(2, 1) + E(2, 2)| < 2, (20) 

\E(l, 1) - E(l, 2) + E{2, 1) + E{2, 2)| < 2, (21) 

I - £(1, 1) + E{1, 2) + E{2, 1) + E{2, 2)| < 2. (22) 
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FIG. 3: (Color online) Plots showing the triangular regions containing values of a and j3 which satisfy inequality (^) (i.e., 
violating the Bell-CHSH inequality). From left to right, the plots are for values of x = —0.8, —0.4, 0, 0.4, 0.8. 



It is not difficult to see that for any given measurement 
directions, at most one of these inequalities can be vio- 
lated. To see this, suppose otherwise and one will find 
that an inequality like the following 



£(1,1) + £(1,2) > 2, 



(23) 



has to hold true. However, this is logically impossible 
as the absolute value of each of the correlation function 
|-^( s i) s 2)| is upper bounded by 1. 

Moreover, by symmetry, the probability of finding 
measurement directions that violate any of these equiv- 
alent inequalities is identical. Hence, the probability of 
finding measurement directions that violate any of the 
Bell-CHSH inequality, with the help of relabeling of mea- 
surement settings and/or outcomes is: 



4 x 



28.3185%. 



(24) 



Identifying nonclassical correlation via linear 
programming 

Linear programming is a convex optimization prob- 
lem H^] which can be efficiently solved on a computer. 
For this kind of optimization problem, both the objec- 
tive function and the constraints are linear in the vari- 
able xeRf In its standard form, a linear program reads 
as M: 



minimize c x, 
subject to Ax = b, 
x >- 0, 



(25a) 
(25b) 
(25c) 



where c S R d , b S R d , A is a d' x d matrix, y represents 
a component- wise inequality and is a d x 1 null vector. 

For the experimental scenarios described in the main 
text, the fc-th party is allowed to perform measure- 
ments in two different measurement directions fij^ , with 
Sfc = 1,2. Let us denote by of} the measurement outcome 
corresponding to the measurement direction fi^ . Now, 
recall that a given correlation admits a locally causal/ 
realistic description if and only if all local measurement 



outcomes are (probabilistically) determined in ad- 
vance solely by the choice of the local setting Sk- Equiva- 
lently, a locally causal correlation is one where the exper- 
imental statistics admit a joint probability distribution 



i 1 



M Jk] 



for all the 2n variables o$ , and which recovers all the 
observed experimental statistics as marg inals [El. For 
example, the joint probability p± n of observing op = 1 
and 02"' = — 1, must be reproducible via: 



Pl „( W = l, N = -i) 
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k=2 
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in n w. 
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±1 



likewise for all the correlation functions, such as 



E(si,s n ) 



E 

k—l,n , 



e n 



\i\ 

°s Pit 



(<M)> 



which must also be reproducible from appropriate linear 
combination of the entries of the joint probability distri- 
bution p\... n (see Ref. p3[ for a more elaborate exposition 
of these constraints). 

Cl early , these constraints can be put in the form of 
Eq. (|25b[) , with x £ R 2 being the entries of pi... n sought 
for, with all (but one) entries of b £ R 3 being the mea- 
sured correlation functions £(si, . . . , Sfe), and with en- 
tries of A giving the coefficients required to reproduce 
each E(si, . . . , s&) from the entries of pi... n . On top of 
the 3" — 1 correlation functions, b must also contain an 
additional entry that is set to 1 ; the same goes for the cor- 
responding row of entries of A, which must all be set to 1. 
This last requirement ensures that the joint probability 
distribution px... n is normalized while the non-negativity 
of the probabilities is taken care of by Eq. (25c). 

The problem of determining if a given correlation ad- 
mits a locally causal description can then be cast as the 
following linear programming feasibility problem: 



minimize 
subject to 



0, 

Ax = b, 

x y 0. 
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For an alternative formulation of the problem as a linear 
programming optimization problem, see Rcf. [p0[. 




0.25 



FIG. 4: (Color online) Probability of finding the measure- 
ment statistics on a pure two-qubit state |$) = cos 6*10)1 10)2 + 
sin 0| 1) 1 1 1)2 to violate the class of Bell-CHSH inequalities. 
Here, the largest Schmidt coefficient cos# quantifies the en- 
tanglement. The values 9 = and 7r/4 give, respectively, the 
pure product state and the maximally entangled two-qubit 
state. Inset: a magnified plot of the crossover region. 



Probability of violation vs entanglement 

For n = 2, a pure quantum state can always be written 
as \ty) = cos#|0)i|0)2 +sin0|l)i|l)2 for some local bases. 
We numerically compute the probability of violation as a 
function of entanglement, as shown in Fig. [|. Note that, 
while the probability of violation is greater with ROM 
in most instances, it also decreases more rapidly to zero 
with decreasing 6. In fact, at 6 « 0.068 ir, the probability 
of violation of |*) is already close to 10~ 5 % with ROM 
whereas the chance is still about 0.26% with RIM (inset, 
Fig. 0). 
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